Three-level mixing and dark states in transport through quantum dots 
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We consider theoretically the transport through the double quantum dot structure of the recent 
experiment of C. Payette et al. [Phys. Rev. Lett. 102, 026808 (2009)] and calculate stationary 
current and shotnoise. Three-level mixing gives rise to a pronounced current suppression effect, the 
character of which charges markedly with bias direction. We discuss these results in connexion with 
the dark states of coherent population trapping in quantum dots. 
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In a recent experiment Payette and co-workers 
studied the transport through a double quantum dot 
(DQD) in which the source-side QD (QDl) had a sin- 
gle electronic level within the transport window, whilst 
the drain-side dot (QD2) possessed three (see Fig. [TJ. 
Gate voltages enabled the position of the former "s-level" 
to be adjusted and thus used as a probe of the second 
QD. Due to non-ellipticity, the levels of QD2 were found 
not to be the familiar Fock-Darwin (FD) levels 0], but 
rather mixtures of them. This gave rise to a distinctive 
feature in the tunneling magnetospectrum consisting of 
an avoided crossing with a central line running through 
it. Strikingly, this central current line was not continu- 
ous as a function of magnetic field, as one might expect, 
but rather showed a strong suppression near the centre 
of the avoided crossing. The authors of Ref. [l[ sug- 
gested a connection between this phenomenon and that 
of the all-electronic coherent population trapping (CPT) 
of Refs. 0, [USUI- It is the aim of this paper to explore 
this connexion further. 

We use a master equation treatment and calculate sta- 
tionary current and shotnoise. We consider a source- 
drain bias direction both as in Ref. [l[ (forward bias), as 
well in the opposite direction (reverse bias). Both bias 
directions yield a current suppression, but as our calcula- 
tions here reveal, the character is rather different in each 
case. In forward bias, the current suppression valley is 
wide (proportional to the mixing energy between the lev- 
els) as observed in the experiment of Ref. [l[ and the shot- 
noise is subPoissonian. In the reverse bias configuration, 
the current suppression valley is narrow (proportional to 
the coupling rate with the leads) and the current statis- 
tics are strongly superPoissonian. We argue that only in 
the latter case does the current blocking mechanism bear 
strong resemblance to coherent population trapping. 



FIG. 1: Double quantum dot with a bias window that in- 
cludes the single probe s-level in QDl and three levels of 
QD2. The depicted bias configuration is as in Ref. [IJ, which 
we describe here as forward bias. In the sequential tunneling 
regime, electron tunneling is described by the rates Tl from 
left lead to QDl, Tr from QD2 to the right lead, and by -y a ; 
a = 0, ± between the dots. 



The Hamiltonian of the first dot reads Hi = e s \s)(s\ with 
\s) denoting the single QDl s-type orbital. Denoting the 
bare FD levels in the second dot as i = 1,2,3, we 
take the Hamiltonian of the second dot to be of the form 



H> 



Eb{\1)(1\ 
+T(|1)(2| 



|3)(3|) 
|2)(1|- 



|2)(3| + |3)(2|), (2) 



where energy Eb = Cb(B — Bq) with cb a constant, 
describes the magnetic field dependence of the FD levels 
(here assumed linear for levels 1 and 3 and constant for 
level 2), and T is the coupling strength between the levels. 
We assume that the coupling between the two dots can 
be described by 



I. MODEL 

We assume strong Coulomb blockade such that at most 
one excess electron can occupy the DQD at any one time 
and write the Hamiltonian of the complete system as 



H — Hi + H2 + H11 



H 



loads 



V. 



(l) 



Hi 2 = nj2(m\ + \i)(s\ 



(3) 



with common coupling parameter f2. Finally, Hamilto- 
nian i?icads = 12k x £kxC lx Ckx describes two standard 
fcrmionic reservoirs (A" = L,R: left, right), and 

v = E (VkLct L \0)(s\ + Y,V kR ct R \0)(i\\ +h.c, (4) 

k \ i=l / 
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with |0) the empty DQD state, describes the coupling of 
the dots to the leads. Note that here we have chosen the 
simplest configuration of parameters — our aim is the 
qualitative understanding of such systems and not the 
quantitative reproduction of the results of Ref. pj . 

The eigenstates of Eq. ([2]) play the determining role in 
the transport through the system; we shall denote them 
|^o) arK1 I^±)j corresponding to cigcncncrgies eg = 



e Q |* Q ); 



and e± = ±y/E% + 2T 2 such that H 2 \^ a ) 
a = 0, ±. The most important of these three states is 
that belonging to eigenvalue zero: 



[*o> 



|*o(£s)> 



(-r|i) + ^ g |2) +r|3» 



(5) 



II. MASTER EQUATION 

All three barriers of the double quantum dot are high, 
and thus a treatment in terms of Fermi's golden rule is 
appropriate for all tunneling in the system. Furthermore, 
we assume large bias such that, for a given bias direction, 
tunneling to/from the leads is unidirectional, and all rel- 
evant Fermi functions are either zero or one. We describe 
tunneling to/from left and right leads with the rates Tl 
and Tr respectively and define T = Tl + Tr. From the 
form of Eq. (|4|) , the right-lead rates are the same for all 
three FD orbitals, and thus also for all three eigenstates 

l*a>. 

Tunnelling between the dots is governed by the matrix 
elements of if 12 and therefore by the overlaps between 



the states in the two dots. Denoting the overlaps of the 
s-level with the QD2 FD states as (s\i) = s,; i = 1,2,3, 
we have, for example, 



1 



(T(s 3 - si) + E B s 2 ) 



(G) 



This eigenstate-overlap clearly vanishes for Eb = T(si — 
S3)/ 's2- In principle, overlaps s; must be determined from 
calculation with orbital wave functions. However, here 
we make the simple assumption that all s» are the same. 
This is justified because the essential feature that (s|^o) 
vanishes remains regardless of the particular values of 
Si. We then set these overlaps to unity, since they can 
be subsumed into the rate 7, to be defined below. The 
squares of the relevant matrix elements arc then 



l(s|*o)| 2 



(si*. 



El 



El + 2T 2 



= 1 



T 2 



± 



2T 



£| + 2T 2 ^JE% + 2T 2 ' 



(7) 

(8) 



Following 0], we then take the hopping rates between 
states s and a to be 



7a =j\(s\H> a }\ 2 L(\e s -E a \,T); a = 0,±, 



(9) 



where we have assumed a Lorentzian broadening of the 
levels, L(x,w) = (1 + (2x/w)^)^ 1 , and 7 = 7(0) sets the 
overall scale for these rates 



With forward bias (as depicted in Fig. [T]) , the Liouvillian (rate matrix) of the system in a basis of populations of 
states ('empty', s, ^q, ,^f+) reads 



£fwd(x) 
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Tr 
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-7+ - 



\ 



(10) 



In reverse bias, the right chemical potential lies above all three levels in QD2, with that on the left lying below the 
QD1 s-level. The Liouvillian for this situation is 



A-ev(x) = 



/ -3T R T L e lx 

-T L - 70 - 7- - 7+ 7o 7- 7+ 

t_r 70 -70 

T R 7- -7_ 

\ T R 7+ -7+ 



\ 



(11) 



Here we have added counting field \ 1° facilitate the cal- 
culation of the current and shotnoise HEH. The density 
matrix itself evolves under the action of the x = Li- 
ouvillian. For example, in the forward bias case we have 

p(t) = rfwd(o)p(t). 



A. Current statistics formalism 

The current statistics of our model can straightfor- 
wardly be calculated using the jump-super-operator for- 
malism of full counting statistics [ll|, 03 ■ The sta- 



3 




tionary density matrix of the system, written as vector 
Iftrtat)) is defined by £(0)|p stat )) = 0. The correspond- 
ing left "Trace vector" is ((cf>o\ = (1,1,1,1,1) such that 
((0o|Pstat)) = 1- Let us define the stationary state "ex- 
pectation value" ((...)) = ((0o | . . . |pstat))) jump super- 



operator J 



£(x)l 



and the pseudo-inverse 



propagator K(z) = Q[z- £(0)] _1 Q with Q = 1 - V 
and V = |/O s tat))((0o|- ln this notation, the current and 
zero-frequency shotnoise read 

{I) = ((J)) (12) 
S = ((J))+2((JK(0)J)). (13) 

We further define the Fano factor as F = S/ (I) . 



III. RESULTS 

In all the following, we set Tl = Tr for ease of pre- 
sentation. Furthermore, in the experiment tunneling 
rates were much smaller that the level-mixing strength, 
T 3> Tr, 7, we will use this fact in various approximate 
results. 
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Figure [5] shows the current through the system as a 
function of magnetic energy Eg and QD1 level position, 
e s . The general structure of the measurements of Ref [l[ 
— an avoided crossing with a line through the middle — 
is reproduced, with current suppression near e s = Eb = 
clearly present. Near this point, the s-level is close to 
resonance with the QD2 state |\&o) an( i, if we ignore con- 
tributions from the other two levels, the current through 
the system may be approximated as 



(I) 



fwd 



7o r J?. 
T_r + 37o ' 



(14) 



The rate 70 is proportional to the matrix clement 
I (s|^o) | 2 which, from Eq. ([7]) is seen to vanish at Eb = 0. 
Within this approximation, the stationary state of the 
system at e s = Eb = is p s t a t = l s )( s li with an electron 
trapped in the s-level due to the vanishing of the ma- 
trix clement. In this approximation, the current at this 
point is zero. From Fig. ^jp, however, it is clear that the 
current is not completely suppressed at e s = Eb = 0, 
but is finite due to the conduction through the other two 
states |^±). This residual current can be estimated as 
ifwd ~ 37r|./(2T 2 ), which need not be negligible. 

The width of current suppression feature at e s = 
can be approximated as follows. Close to Eb = 0, 70 
is small, and Eq. (TTJJ can be further approximated as 
Iiwd ~ 7o- On the other hand, far from Eb = 0, the 
current saturates to the constant value if w d ~ 1^7/(37 + 
r^). The value of the magnetic field at the point where 
these two behaviours cross can be found by setting the 



FIG. 2: Current through the DQD in the forward bias con- 
figuration: (a) density plot as a function of magnetic energy 
Eb and s-level position e a . (b) sections at e s = (black), 
e s = — 25Tr (blue), and e s = 25Tr (green). Also plotted 
(dashed) is the approximate current of Eq. (|14|l (for e 3 — 0) . 
A strong current suppression is observed around the point 
e s = Eb = 0, but note that a small current does flow at this 
point, however. Further parameters are: 7 = Fl — Tr and 
T= WTr. 



two limiting values equal, and solving for Eb- Equating 
this value to half the width of the current suppression 
valley we find 



Wlwd = J —T, 



(15) 



which shows the width of the current suppression valley 
to be proportional to the level-mixing energy T. 

The shotnoise Fano factor for this bias direction is 
shown in Fig. [3] Especially evident is that the Fano factor 
is everywhere less then (or equal to) unity, corresponding 
to the familiar subPoissonian statistics of anti-bunched 
electron transfer. Again assuming that only the central 
resonance determines the transport in the neighbourhood 
of the current suppression, we can approximate 



F, 



fwd 



F 2 

1 R 



2r i?7o + 5 7o 2 



(Tr + 3 7o ) 2 



(16) 



which is clearly always less than or equal to unity. Along 
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FIG. 3: As Fig. [2j but here the shotnoise Fano factor is 
shown. Along the resonant lines of high current, the shotnoise 
is subPoissonian with a Fano factor F ~ 1/2. Away from 
resonance, as well as around the central current suppression 
region, the shotnoise approaches the Poissonian value F = 1 
from below. 



FIG. 4: Current in the reverse bias, parameters as Fig. [2] 
Here, the level structure of the second dot is not resolved. 
Rather, two large current peaks are observed. Along the 
Eb = axis, the current is completely suppressed, irrespec- 
tive of the s-level position. This is attributed to the formation 
of the CPT dark state of Eq. (fT7)l . 



the central resonance (e s = 0), the Fano factor reaches 
a maximum value F f=a 1 — 6jTr/T 2 at Eb = and a 
limiting value of F w {5j 2 + 2jT R + T 2 R )/(T R + 3j) 2 for 
large Eb along the e s = line. 



B. Reverse bias 

Figure 2] shows the current with the source-drain bias 
in the opposite direction. Once again, the current shows a 
suppression at Eb — 0, but unlike the forward bias case, 
this suppression extends for all positions of the s-level. 
The second significant feature of this suppression is that 
the current is exactly zero at Eb = 0, even for T/7 finite. 
It is easily shown that for Eb = 0, the stationary density 
matrix of the system is p sta t = | , I'o(0))(\E f o(0)|, which 
clearly shows that in the long time limit, an electron is 
trapped in the DQD in the pure state 



l*o(0)) 



1 



V2 



(|3)-|1» 



(17) 



We can obtain approximate expressions for the cur- 
rent as follows. If Tr and 7 are of the same order-of- 



magnitude, then, for any choice of e s , at least one of 7 Q 
will be much smaller than Tr. We can then write the 
current as 



I* 



3707+7- 



7o7+ + 7o7- + 7+7- 



(18) 



Concentrating about the point e s — 0, we can say that 
near Eb = there is always a regime in which 70 is 
the smallest rate (matrix element disappears). In this 
case the current becomes / rov ~ 37o- This expression 
describes the form on the sharp dip in the current about 
Eb = 0. Further away from Eb = 0, 70 becomes the 
largest of the 7, rates since its matrix element returns 
to a typical non-suppressed value, and the other two 
rates are off-resonant. In this case, the current becomes 
J rcv w 37+7-/(7+ +7-). An estimate of the width of the 
current-suppression feature can then be obtained from 
the cross over between these two behaviours which oc- 
curs when 70 = 7+ (NB: for e s = 0, 7_ < 7+). Solving 
for Eb, we find the width to be 



•mv 



6 + 4V2T R , 



(19) 
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FIG. 5: Shotnoise Fano factor in the reverse bias, parameters 
as Fig. [3] In contrast to the forward bias case, the shotnoise 
is almost everywhere superPoissonian. For Eb = 0, the Fano 
factor is exactly 5. 



which is proportional to the broadening induced by the 
contacts. 

The corresponding Fano factor is shown in Fig. [5] In 
strong contrast to the forward bias case, the noise here is 
almost everywhere superPoissonian, and in particular in 
the neighbourhood of the current suppression. Analytic 
expressions for F are, in this case unwieldy. However, 
without any further approximation, we find that at Eb = 
the Fano factor is simply F = 5 independent of all 
further parameters. 



For forward bias, near the point es = Eb = 0, con- 
duction comes through three channels which are weakly 
transmitting: two (±) on account of their distance in en- 
ergy from resonance with the probe level at e s , and one 
(0) on account of the vanishing of the matrix clement for 
hopping between the two dots. In this case, the steady 
state of the DQD is approximately that of an electron 
trapped in the QD1 probe level p s tat ~ \s) (s\. This trap- 
ping is not exact, however, and a current still flows at 
«s = Eb = due to conduction through the off-resonant 
channels. As a natural consequence of having a set of 
weakly transmitting channels, the statistics are subPois- 
sonian. This situation resembles somewhat the isospin 
blockade of Ref. [T3]. 

In contrast, the current in the reverse direction for 
Eb = is exactly zero — not just for e a = 0, but ir- 
respective of probe-level position. In this case, the dot 
electron is trapped in the state: 2 -1 / 2 (|3) — |1)). This 
is a pure superposition state and is directly analogous to 
the dark state of the triple quantum dot CPT Q. As in 
the triple dot case, the corresponding current statistics 
are superPoissonian. This may be understood in terms 
of the dynamical channel blockade [l4j], since we have 
one weakly transmitting channel (that associated with 
the dark state), and two normally conducting ones. 

A further distinguishing feature between these two 
blockade situations is that the width of the forward bias 
suppression valley is proportional to the mixing ampli- 
tude T (the large energy scale in the model) , whereas that 
in reverse bias is proportional to the lead-coupling rate 
r# (the small energy scale). We also mention that the 
reverse-bias suppression is robust if we increase the inter- 
dot coupling Q, whereas the forward-bias feature washes 
out as the three resonances start to overlap. 

The reverse-bias CPT effect described here should be 
more robust with regards to dephasing induced by e.g. 
background charge fluctuations than the dark state of 
the triple QD since it is substantially more localised. No 
interdot coherence is required. Furthermore, the ease 
in tuning the magnetic field to precisely locate the dark 
state, as compared with the gate voltages in the triple 
QD, makes this an excellent set-up for the further study 
of CPT and dark states in mesoscopic transport. 



IV. DISCUSSION 
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